Single-copy entanglement in a gapped quantum spin chain 
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The single-copy entanglement of a given many-body system is defined [J. Eisert and M. Cramer, 
Phys. Rev. A. 72, 042f 12 (2005)] as the maximal entanglement deterministically distillable from 
a bipartition of a single specimen of that system. For critical (gapless) spin chains, it was recently 
shown that this is exactly half the von Neumann entropy [R. Oriis, J. I. Latorre, J. Eisert, and 
M. Cramer, Phys. Rev. A 73, 060303(R) (2006)], itself defined as the entanglement distillable in 
the asymptotic limit — i. e. given an infinite number of copies of the system. It is an open question 
as to what the equivalent behaviour for gapped systems is. In this paper, I show that for the 
paradigmatic spin-S AfHeck-Kennedy-Lieb-Tasaki chain (the archetypal gapped chain), the single- 
copy entanglement is equal to the von Neumann entropy: i.e. all the entanglement present may be 
distilled from a single specimen. 

PACS numbers: 75.10.Pq, 03.65.Ud, 03.67.-a 



The amount of entanglement naturally present in the 
ground states of many-body quantum systems is of great 
interest to both the quantum information and condensed 
matter communities [3]. For bipartite pure states the 
unique asymptotically continuous entanglement mono- 
tone is the von Neumann entropy the amount of 
entanglement present in the ground states of spin chains 
and other many-body sytems is often quantified by calcu- 
lating this quantity with respect to a particular biparti- 
tion. It is asymptotically equal to the ratio of the number 
of EPR singlets one may distill from a given system to 
the number of copies of that system — that is, given an 
infinite number of copies of the system Of course 
in realistic, physical situations one might only have ac- 
cess to a single specimen of a particular system, and the 
von Neumann entropy in this case gives an upper bound 
to the distillable entanglement. An interesting question 
then arising is: how much entanglement can be deter- 
ministically distilled from a single copy of a system? 

This was studied originally by Lo and Popescu , who 
considered various optimal strategies for the interconver- 
sion of entangled states by loj;al operations. It was sub- 
sequently shown by Nielsen J] that deterministic trans- 
formation of one state 1-0) to another \(f>) using local op- 
erations and classical communication is possible if the 
nonincreasingly ordered Schmidt coefficients of the final 
state majorize those of the initial state (A^ -< A^); this 
was then generalised to a strategy for converting arbi- 
trary bipartite states by Vidal [8, 9]. 

More recently, Eisert and Cramer [1] and Orus et al. [2] 
defined the single-copy entanglement as the number of 
singlets deterministically distillable from a single speci- 
men of a given system, and showed that for the case of 
gapless quantum spin chains close to criticality this is 
exactly half the von Neumann entropy: that is, half the 
entanglement present may he distilled in a single process 
[4l| . However, a large class of spin chains of interest are 
gapped (i.e. there is an energy gap between the ground 



and first excited states), which have substantially dif- 
ferent entanglement properties: in the gapless case, the 
block entropy (the entanglement of a block of contiguous 
spins with the remainder) is known to be logarithmically 
divergent in the block l eng th, and in the gapped case, 
it saturates 13, U , [13, [3, [31 • It is this short-ranged 
nature of entanglement i n g apped systems which gives 

MM 



rise to area laws 15, 16 



It would therefore be 
interesting to know whether gapped chains also exhibit 
different behaviour in the context of single specimens. 
While this remains an open problem in general, I shall 
in this paper provide an example giving indications that 
this is indeed the case. Specifically, I shall show that the 
celebrated spin- 5 Affieck-Kennedy-Lieb-Tasaki (AKLT) 
chain [3, 2y, 2lj has a single-copy entanglement equal to 
the von Neumann entropy for all S: i.e. all the entangle- 
ment present in the spin chain may be distilled in a single 
process. This spin chain is a paradigm in condensed mat- 
ter physics, and was the first evidence of the veracity of 
the Haldane conjecture [H, HI]. 

Recent years have seen a resurgence of interest in this 
state, partly due to its role as the simplest of the ma- 
trix product states, which in general have been shown 



to efficiently simulate many ID systems [2J, |25|, |26!, |27| 
and may be used as a variational set in density matrix 
renormalisation group calculations [13, H^] ■ It has also 
attracted interest in the context of quantum informa- 
tion theory, due to its interesting entanglement proper- 
ties (s^, [m [sl], and the fact that all stabilizer states 
(including cluster states) have an interpretation in terms 
of a VBS state ^33,,i34,]- 

Single-copy entanglement. — The single-copy entangle- 
ment (with respect to a particular bipartition) is defined 
[ij as the maximal number of singlets one can determinis- 
tically distill from a single-copy of a specimen in a single 
process; i.e. the single-copy entanglement is Ei = logM, 
if M is the largest m for which the transformation 
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is possible under local operations and classical com- 
munication with unit probability (where |V'm) = 
l/i/m^^j^ I*,*), the maximally entangled state of di- 
mension m X m). This is the equivalent number of sin- 
glets ^2|. Applying Nielsen's majorization criterion to a 
bipartition of a spin-S* chain into a block of length L and 
the remainder, this holds if and only if 
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5]«i<^foralli^e[l,M], 



k=l 



M 



(2) 



where {aj, . . . ,a 



(2S+1) 



, } arc the nonincreasing eigenval- 



ues of the reduced density matrix of the block of length 
L. As pointed out by Eisert and Cramer [H, this is obvi- 
ously equivalent to the criterion that a\ < 1/M (ie. the 
majorization reduction must be more mixed than the 
maximally entangled state). One can then define the 
single-copy entanglement 



£;i(p)=log l/ai 



logaj. 



(3) 



Critical, gapless systems. — Using the machinery of 
conformal field theory, Oriis et al. found that for all trans- 
lationally invariant quantum spin systems that can be 
mapped onto an isotropic, quadratic system of fermions 
(via the Jordan- Wigner transformation), the single-copy 
entanglement of a block of length L is exactly half the 
von Neumann entropy in the thermodynamic (L —^ oo) 
limit: both logarithmically diverge with L. Explicity, 

E,{p,) = ^lnL-^^ + 0{l/L), (4) 

where c is the conformal field theoretic central charge. 

Noncritical, gapped systems. — It is known that for 
gapped systerns, the block entropy saturates to a con- 
stant bound [lO, 0, [sSl- This qualitatively different 
entropic behaviour would suggest that the single-copy 
entanglement might also have substantially different be- 
haviour. As a first step towards understanding this, it 
would be instructive to consider the AKLT chain, the 
first example of a chain satisfying the Haldane conjecture 
igl . 20l . 21| . This chain consists of spins of magnitude 
S, with two spin-S'/2's at either end. In the Schwinger 
boson representation this is written [35 1 
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|VBS) = Yliajbl, - bjal.f |0) 



(5) 



i=0 



where aj, bj are bosonic operators, and the spin operators 
are defined as 5+ = ajb^, 5"^^ = OifoJ, Sf = (a|ai-6|5i)/2, 
with the constraint a|a,i + bjbi — 2S. For general S, this 
is the unique ground state of the Hamiltonian 

N-l 2S 

H=Y1 ^J^L+i + '"O.l + T^NM+l, (6) 

3 = 1 J=S+1 



o — do- 



On — oo 



FIG. 1: The VBS state has an interpretation in terms of 
'bonds' between spins, where each bond is a singlet of two 
spin-S'/2's, and at each bulk spin the state is projected (dotted 
circles) to the symmetric space of total spin. The single-copy 
entanglement measured here is with respect to the partition- 
ing into a block of length L and the remainder. 



where the operator Pjj+i projects bond spins j, j+1 onto 
the (symmetric) subspace of total spin J, and the Aj are 
arbitrary positive coefficients. The boundary terms ttq,!, 
T^N.N+i are similarly defined to project the end spin S 
and S/2 onto the total spin J subspace. The isotropic 
case with no boundary operators has also been studied for 
the case = 1 [sj , and the edge effects decay exponen- 
tially; similarly the effect of studying periodic boundary 
conditions in Eq. ([6]) as opposed to open boundary con- 
ditions also decays exponentially ^ . This state may be 
interpreted in terms of 'bonds' between spins (Fig. [1]), 
with each bond being a singlet, and is often referred to 
as a valence bond solid. 

Coherent spin state approach to the VBS. — An alter- 
native representation of state (O is given in terms of 
coherent spin states. The well-known oscillator coherent 
state is defined to be an eigenstate of the annihilation 
operator a; by analogy, the coherent spin state is defined 
as an eigenstate of the spin raising operator 5+ [s^] 



I 



(l + lA^P) 



2\2S 



exp(/iS' ) 15, S") , 



(7) 



where \S,m) is the spin state with (S^) = S{S + 1) and 
{Sz) = 1TI- Parametrizing this with /i = e"^ tan(0/2) gives 



E 




S, m) , (8) 



where {u,v) := {e"!'/^ cos{e /2),e-'*/^ sm{e /2)). This 
state has a clear geometric interpretation: the state 
|0, (j)) may be represented by the unit vector Q, = 
(OjCf)) {i.e. a, point on the unit sphere). Therefore 
the overlap between two such states may be found 
geometrically as (6*, 0161', 0') = [cos(6i/2) cos(6l72) -|- 
sin(6'/2) sin(6'72) exp(i(0 - 0'))]^^, and thus 



1 + vt-n' 



(9) 



Comparing Eq. ([5]) to the Schwinger boson representa- 
tion of the S^, Sz eigenstates 

\S.m)^ ;7 \l 10), (10) 

^{S -\- my.yS — my. 
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|2S 



it is clear that one can obtain an alternate representation 
of a state by making the replacement — > u, 6^ — s- v, 
a d/du, b d/dv, and multiplying by ^ (25)! when 
operators occur in pairs. 

Therefore, denoting the coherent state \n) one may 
write [3i,[3i (r2|VBS) : 
and hence | (17|VBS) p 

n^ij) i2S)\{l~'flk-flk+i)^/2^- This approach was used in 
Ref. ^Mi] to calculate all two-spin correlation functions, 
and more recently in Ref. [32| to calculate the block 
entropy. 

Single-copy entanglement of the VBS. — Using the co- 
herent spin state approach, one may calculate the den- 
sity matrix oi L contiguous spins. The trace of 
any operator in this representation is tr ^ = (25* -|- 
l)/47r / dfi (ri|y^|ri) and therefore one obtains 



PL = 



n 

m 



dn (l]j|VBS) (VBSIOj) 



47r (VBS I VBS) 



i=l 4tt 



Ylk=l TkM+l 



r t-tL dilj t-tL — 1 ry-i ' 

J llj = l —[lk=l ^k,k+l 

where the transfer matrix Tk,k+i (1 ^ ^k ■ ^k+i)^ /"^^ 
and in the first line I have omitted numerical factors. It 
was found in Ref. [s^l that this is independent of the 
length of the total chain N, and therefore without loss of 
generality one can set L = N. Following the methods of 
Refs. [s^, [IHl , one may find the eigenvalues of this matrix 
(replacing S S/2 for the end spins) 

The following decomposition in terms of Legendre 
polynomials may be used 



l + x 



SlSl 



1=0 



[s ~ ly.is + 1 + ly. 



from whence it follows that 

PL 



J2x{l)L-^Ii{so-SL+i), (12) 



where A(0 {-Y Sl{S+iy/{S -iy.{S + 1 + iy und Ii{X) 
is an Ith order polynomial in X determined recursively 
through 



2j+3 



AX 



(5 + J + 2)2 \j + 
J 2j + 3 fS- j + 1 



j + 12j-l \S + j + 2 



I,-iiX) (13) 



with Io{X) = l/47r, Ii{X) = 3X/47r(S'/2+ 1)^. These 
polynomials form a complete set of isotropic, two-site 
tensor operators. 

The eigenvalues of the density matrix may be found 
using this method, and summed using the standard for- 
mula for the von Neumann entropy [S{p) — — tr p log p\ to 



find the block entropy for general S', which was found to 
approach 2 log(5'-|-l) exponentially fast in L (the thermo- 
dynamic limit) (3^ . confirming the conjecture by Vidal et 
al. that the block entropy of a gapped integer spin chain 
reaches saturation for all S . For the purposes of this 
paper, only the largest eigenvalue is required. 

The density matrix is diagonal in the basis of the total 
spin of spins and L+1. These spins, of course, add up 
to several multiplets in the usual manner of spin addition, 
and so there will be degeneracy in the eigenvalues. One 
thus requires the largest value of 



(Fa) =tr {PaPL}, 



(14) 



where is the projector on to the subspace of total 
spin a. This multiplet distribution (i.e. the eigenvalues 
multiplied by their weight) is given by [s^ 



(P,) = 47r(2a - 



l)V^^±^±^I^A(,)-^/,[X(.)], 



(5 + 1)1(5 + 1)! 



(15) 



where X(o-) := Sq-sl+i = cr(cr + l)/2-5/2(5/2 + l). The 
calculation of the single-copy entanglement only requires 
the eigenvalues, and thus one omits the weights {2a + 1). 

This distribution is found recursively, determined by 
the coefficients /j[X(cr)]. The largest values are given by 
the case cr = 5 (ct = 0) for L even (odd). For even L, the 
required value is 



{Ps)^i2S + l)J2 



2i + i 



L + l 



and thus the largest eigenvalue is 

2j + l 

(5- 



(16) 



(17) 



A-E#^A(,)- 

j=Q 

1^|i + EA(j)^+1(2, + 1)|, (18) 



is 

which gives the single-copy entanglement: 
El = -logAi 



= 21og(5+l)-log<^ 1 + EA(j)'^+'(2j + 1) \ (19) 



= 2 iog(5 + 1) - log A(o)^+i 1 1 + E ^Ijirr (2j + 1) 

(20) 

Since A(j) > A(j + 1) for all j, and A(0) — 1, it is clear 
that in the thermodynamic case L — > 00 (as considered 
in the critical, gapless case 0), this becomes 



£;i^21og(5+l), 



(21) 
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which is exactly equal to the von Neumann entropy, as 
found by Katsura et al. [s^l ■ The proof for L odd follows 
analogously. It therefore is the case that all the entan- 
glement present in the VBS state (the ground state of 
the gapped spin-S" AKLT Hamiltonian) may be distilled 
from a single copy: one can distill with certainty a max- 
imally entangled state, the dimension of which is related 
to S. This would appear to have an intuitive explana- 
tion in terms of the valence bond picture of the state: 
the entanglement between a block and the remainder of 
the chain is related to the number of bonds 'cut' by the 
boundary (indeed, this is similar to the reasoning behind 
area laws [isl \wi 17 . although the analogy is not 
strict in this case, since the entanglements of formation 
and distillation are only equal in the asymptotic limit), 
and is further evidence of the qualitatively different be- 
haviour of gapped chains to gapless chains. One should 
contrast this with the critical case 0, from whence one 
can distill with certainty (in the L — > oo limit) a maxi- 
mally entangled state of arbitrary dimension {i.e. an in- 
finite single-copy entanglement); the crucial, qualitative 
difference is that this is still only half the total amount 
of entanglement present. 

Summary. — It has been demonstrated that all the 
entanglement present in the VBS ground state of the 
gapped AKLT Hamiltonian for arbitrary S may be dis- 
tilled with certainty in a single process. This quali- 
tative difference from the behaviour of gapless, critical 
chains provides evidence that the entanglement present 
in gapped systems is of a fundamentally different nature. 
An open problem is whether the single-copy entangle- 
ment saturates to the von Neuman entropy for all gapped 
quantum systems. 
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